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EXECUTIVE  SUMMARY 


As  the  title  suggests  this  document  provides  brief  instructions  for  a  set 
of  computer  subroutines  called  SIMBAT.  SIMBAT  can  unconditionally  and  condi¬ 
tionally  simulate  random  directional  ocean  wave  properties.  Given  an  estimate 
of  the  directional  wave  spectrum,  the  program  calculates  elevations,  kinematics, 
and  pressures  in  the  random  wave  field.  These  calculations  can  be  conditioned 
on  particular  kinematic  or  surface  elevation  time  series. 

Conditional  simulation  of  wave  property  time  series  statistically  consis¬ 
tent  with  a  specified  measurement  set,  provides  a  very  powerful  approach  to 
certain  ocean  engineering  problems.  The  usual  computer  simulation  of  waves 
satisfying  a  specified  model  for  the  directional  spectral  density  suffers  from 
a  serious  practical  defect  if  one  is  primarily  interested  in  producing  very 
large  waves.  Most  simulations  produce  only  average  waves  unless  the  simulation 
is  run  for  a  very,  very  long  time.  SIMBAT  allows  for  the  inclusion  of  a  large 
wave  profile  or  wave  group  to  be  embedded  into  the  wave  train,  resulting  in 
very  short  computer  simulations. 

The  report  describes  basic  wave  properties  in  their  complex  form,  describes 
the  program  SIMBAT,  and  explains  in  detail  the  development  of  the  Legendre 
polynomials  for  storage  of  the  large  volume  of  wave  kinematic  data  generated. 

This  contract  report  was  prepared  by  Dr.  Leon  Borgman,  professor  of 
Statistics  and  Geology  at  the  University  of  Wyoming,  working  for  the  Naval 
Civil  Engineering  Laboratory  through  his  statistical  consulting  firm,  Leon  E. 
Borgman,  Inc.  The  work  was  principally  funded  by  the  Mineral  Management 
Service  through  Charles  Smith  of  the  Technology  Assessment  &  Research  Branch. 
Additional  work  in  fiscal  year  1990  is  planned  under  Naval  Facilities  Engineer¬ 
ing  Command  funding  for  testing,  modifying,  and  annotating  SIMBAT. 
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1.0  GENERAL  COORDINATE  SYSTEM 

The  ocean  wave  kinematics  will  be  referenced  to  a  general  horizontal 
coordinate  system.  All  horizontal  coordinate  axes  are  established  within 
navigation  headings  measured  clockwise  from  true  north. 

8  =  direction  of  positive  x-axis 

x 

0^  =  direction  of  positive  y-axis  (1) 

|0  -  0  I  *  90° 

1  x  y1 

Let  the  vertical  axis  z  be  zero  at  mean  water  level  and  positive  down¬ 
ward. 

The  direction  of  travel  of  a  wave  is  8  in  navigation  heading.  The 

wave  is  traveling  toward  direction  0  if  (JQ  ®  1  and  is  coming  from 

direction  0  if  0  *  -1. 

o 

2.0  BASIC  WAVE  PROPERTIES 

Eight  wave  properties  are  of  interest.  In  terms  of  real  functions, 
there  are 

(1)  The  water  level  elevation: 

n(x»y.t)  *  a  cos  {0q  k  [x  cos(0  -  0x)  +  y  cos(0  -  0  )] 

-  2i»ft  -  (2) 

(2)  The  components  of  water  perticle  velocity: 

Vx(x,y,z,t) 

Vy(x,y,z,t) 

•  cos  [8^  k{x  cos(0-8x)  +  y  cos(0-0^)) 

-  2irft  -  *]  (3) 


a(2irf)  c-osh^;-^ 
s  mh(kd) 


6  cos(9-8  ) 

o  x 

8  cos(8-8  ) 

o  y 


1 


(3) 


V  (x,y,z,t)  = 

The  components 


,iat{£;;)1  si"'6okix  c°s<9-v 

+  y  cos(8-0  ))  -  2irft  -  $] 
of  water  particle  acceleration: 


ax(x,y,z,t) 
ay(x,y,z,t) 

•  sin[B  k{x  cos(0-0  )  +  y  cos(0-0  )} 

OX  y 

-  2 it f t  -  4] 

az(x,y,z,t )  =  -  a(2irf)2  cos[BQk{x  cos(0-8x) 

+  y  cos(9-0y)j  -  2ifft  -  0] 


B  cos (0-0  ) 
o  x 

B  cos(0-0  ) 

o  y 


(4)  The  water  pressure  anomaly  (plus  and  minus  about  hydrostatic 
pressure) : 

P(x.y.z.t)  =  apg  cos[B0k{x  cos(0-0x) 

+  y  cos(0-0^)|  -  2irft  +  $] 

In  these  formulas 

a  =  wave  amplitude 
f  *  wave  frequency 
d  *  water  depth 

k  »  wave  number  =  2x/wave length 
%  *  wave  phase 
p  =  water  density 
g  =  acceleration  due  to  gravity 


(4) 


(5) 


(6) 


(7) 
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3.0  WAVE  PROPERTIES  IN  COMPLEX  FORM 


Through  the  use  of  the  complex  form  of  cos  a  and  sin  a  where 

cos  a  =  exp(ia)  ^exp(-ia)  (g) 

ain  a  =  -_,exP(-ia)  (9) 

All  of  the  wave  properties  listed  above  can  be  expressed  in  the  form: 
i$ 

B(f)  =  G(z)  T(f)  H(0)  exp[-i0ok[x  cos(0-0x) 

+  y  cos(0-9y)}]  exp(i2irft)  (10) 

for  positive  f.  The  original  real-valued  wave  time  property  equals 
B(f)  +  B( f )*  where 


B*(  f ) 

3 

complex  conjugate  of  B(f) 

(ii) 

The 

functions  1 

3,  T,  and  H  for  each  wave  property  are: 

(1) 

Water 

level  elevation: 

G(z) 

=  T(f)  =  H(0) 

5  1.0 

(12) 

(2) 

Velocity: 

coshfk(d-z) 1 

,  -kz  , 
_  {«  + 

-k(2d-z) . 

e  I 

for  V  and  V 
x  y 

sinh(kd) 

{1  - 

-2kd, 
e  ) 

G(z) 

sinhf k(d-z) ] 

.  -kz 
_  l« 

-k( 2d-z) . 
e  i 

(13) 

for  V 

z 

sinh(kd) 

- 

-2kd. 

■ 

e  J 

'  2irf  , 

for  V  and 

V 

T(f) 

2irfi  , 

. 

X 

for  V 

z 

y 

(14) 

B  cos(8-9  ) 

,  for 

V 

X 

X 

H(0) 

= 

B  cos(0-0  ) 

,  for 

V 

(15) 

y 

y 

1.0 

,  for 

V 

3 


(3)  Acceleration: 
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Also  let  k  be  replaced  by  k  where 

m 

(2vfh  =  g  k  tanh(k  d)  (22) 

m  mm 

The  sequence,  defined  above  for  1  £  m  <  N/2,  can  be  extended  to  N/2  <  m 
£  N-l  by  requiring  that  (in  analogy  to  Eq  10) 

B[ (N-m)Af ]  =  B(m  Af)*  (23) 

where  B(f)  is  the  general  complex-valued  wave  property  defined  previously. 
The  sum  of  these  wave  forms  over  0  £  m  £  N-l  gives  a  discrete  Fourier 
transform  version  of  the  wave  properties.  Here,  it  is  assumed  that 
=  0  for  m  =  0  and  m  =  N'2.  The  m  =  0  value  is  the  mean  or  DC  component. 
Taking  it  as  zero  guarantees  that  the  wave  property  oscillates  about 
zero.  The  value  at  m  -  N/2  is  a  very  high  frequency  component  at  the 
Nyquest  frequency.  The  length  of  the  time  series,  N,  and  the  time 
increment,  At,  can  always  be  selected  so  that  there  is  no  energy  at  f  ,, 

=  (N/2)  Af. 

Then, 


{wave  property] 
at  t  *  n  At  J 


N-l 


I 


m=0 


1 


A  .  G  T 
mj  m  it 


H. 

J 


-iB  k  (x  cos(9,-0  )  +  y  cos(0.-9  )) 
.  e  o  ml  j  x'  J  y 


i2n 

e 


mn 

N 


.  (24) 

This  represents  a  summing  of  many  waves,  each  with  their  own  frequency, 
phase,  and  direction. 

The  last  equation  provides  the  general  procedure  for 
frequency-domain  wave  simulation.  The  quantity  in  the  bracket  is 
computed  for  1  S  m  <  N/2.  Usually  this  is  only  necessary  for  a 
relatively  small  subset  of  the  interval,  say 

0  <  mg  S  m  S  mL  <  N/2  (25) 
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or 


(26) 


NUM  =  mT  -  m_  +  1 

li  D 

frequency  increments.  Then  the  rest  of  the  coefficients  for  0  <,  m  S  N/2 
are  set  to  zero.  The  values  for  N/2  <  m  5  N-l  are  the  complex  conjugate 
of  those  in  the  left  half  of  the  sequence.  That  is, 

{coefficient  at  N-mJ  =  {coefficient  at  mj*  (27) 

PROGRAM  LAYOUT 

The  program  has  seven  options  in  addition  to  exit  and  help  options. 
These  are: 

(1)  Option  No.  1. 

A  complex-valued  matrix  of  wave  amplitudes  in  the  form: 

*/  ...  ,  .  \  i$(m,j) 

A(m,j)  =  p(m,j)  e 

is  simulated  by  frequency-domain  computations.  Here,  p  is  the  wave 
amplitude,  £  is  the  phase.  The  m-index  ranges  over  a  regular  grid  of 
frequencies,  and  the  j-index  ranges  over  direction  of  wave  travel.  This 
option  gives  an  unconditional  simulation. 

(2)  Option  No.  2. 

A  conditional  simulation  is  a  time  series  simulation  which  is 
forced  to  agree  with  a  specified  initial  data  time  series  segment,  while 
maintaining  appropriate  correlated  randomness.  The  program  SIMBAT 
develops  conditional  simulation  by  several  methods,  all  based  on  first 
producing  a  conditional  simulation  of  the  A(m,j)  described  under  the 
first  option.  These  A(m,j)  however  are  conditionally  simulated,  rather 
than  unconditionally  simulated.  The  production  of  such  a  set  of 
conditional  simulation  requires  a  number  of  pre-computed  arrays.  Option 
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No.  2  develops  these  input  arrays.  Thus,  it  is  a  pre-processor  to  the 
various  options  which  subsequently  compute  the  actual  conditional 
simulations . 

(3)  Option  No.  3. 

This  option  uses  the  output  from  the  previous  option,  for  the  case 
where  the  conditioning  interval  is  shorter  than  the  full  time  series,  to 
compute  a  conditional  simulation  of  the  A(m,j)  complex  matrix  of  wave 
amplitudes . 

(4)  Option  No.  4. 

This  option  is  the  same  as  Option  No.  3  above,  except  that  the 
conditioning  interval  is  a  full  time  series.  That  is,  a  measured  time 
series  of  length,  N,  is  used  to  develop  the  A(m,j)  complex-valued  wave 
amplitude,  which  in  turn  may  be  used  to  simulate  in  a  later  option,  time 

series  of  length,  N. 

(5)  Option  No.  5. 

This  option  uses  the  complex  amplitude  table  of  A(m,j)  values  to 
generate  full  time  series  (of  length  N)  for  various  wave  properties  as 
specified.  It  basically  is  designed  to  be  useful  for  the  case  where 
only  a  few  time  series  (say  20  or  less)  are  needed. 

(6)  Option  No.  6. 

This  option  is  the  fastest  way  to  develop  velocities  and  accelera¬ 
tions  at  many  load  points  throughout  a  complex  structure.  It  provides 
orthogonal  polynomial  coefficients  for  a  Legendre  expansion  of  each  of 

eight  wave  properties  (ti,V,V,V,a,a,a,  pressure)  within  a 

x  y  2  x  y  z 

region  (x^Dj  <  x  <  x^+Dj),  (yQ-D2  <  y  <  yQ+D2),  and  0  <  z  <  d.  A  dif¬ 
ferent,  (optional)  set  of  coefficients  are  provided  at  time  step,  as 
stored  on  a  file.  Option  No.  6  computes  these  coefficients  and  stores 
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them  on  a  master  file.  The  user  then  reads  the  master  file,  time-step 
by  time-step,  and  uses  coefficients  at  that  time  step  to  generate  all 
the  wave  kinematics  at  the  various  (x,y,z)  locations  throughout  the  struc¬ 
ture.  Then  the  user  reads  the  next  coefficient  set  at  that  time  step, 
and  so  forth. 

(7)  Option  No.  7. 

This  option  provides  a  short  list  of  wave  amplitudes,  phases, 
frequencies,  wave  numbers,  and  travel  directions  which  give  a  wave  train 
that  approximates  the  full  wave  field  represented  by  the  A(m,j)  matrix. 

It  will  not  exactly  agree  with  the  condition  set,  but  will  be  somewhat 
near.  It  is  probably  slower  also,  in  use,  than  Option  No.  6  above. 

An  overall  flow  chart  is  shown  in  Figure  1.  The  delta  stretching 
is  applied  by  the  user  as  based  on  having  the  sea  surface  elevation  and 
the  wave  kinematics  at  the  same  (x,y)  location.  Other  spectral  models 
and  spreading  functions  will  also  be  coded  into  the  package.  Currently, 
the  Ochl-Hubble  and  Gaussian  spreading  function  are  implemented.  However, 
it  is  easy  to  insert  other  choices  as  subroutines  into  the  structure. 

It  is  anticipated  that  many  further  enhancements  and  modifictions 
will  be  introduced  during  July  and  August  1988  as  a  natural  outgrowth  of 
experience  and  further  needs  as  the  program  is  used  in  applications. 

6.0  SIMBAT  SEPARATE -MODULE  PACKAGE 

The  SIMBAT  simulation  package  is  also  provided  in  separate  modules 
of  relatively  small  size  for  computation  on  larger  microcomputers.  The 

modules  are: 

1.  SIM125  --  options  1,  2,  and  5 

2.  SIM3  --  option  3 

3.  SIM4  --  option  4 

4.  SIM6  --  option  6 

opt ion  7 


5.  SIM7 


Before  running  program:  Set  parameters  and  compile,  either  for  unconditional  simulation 
(option  no.  1)  or  conditional  simulations  (option  no.  2) 


Figure  1.  Overall  flow  chart  for  SIMBAT.FOR. 
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Figure  1.  Overall  flow  chart  for  SIMBAT.FOR  (continued). 


Each  module  Is  used  by  executing  the  module  and  exiting  with  storage  of 
the  results  on  a  user-selected  file.  The  set  of  complex  wave  amplitudes 

{A(m, j)  ;  0  S  m  S  N-l,  0  i  j  i  J-1J 

Each  option  can  be  classified  relative  to  its  relation  to  A(m,j)  as 
shown  in  Figure  2. 


Option  No.  1 


Option  No.  3 


Option  No.  2 
conditional 
simulation 
preprocessor 


<- 


Option  No.  4 


unconditional 

simulation 


conditional 
simulation 
L  <  N 


conditional 
simulation 
L  -  N 


produce  a  finite 
number  of  time  series 


Option  No.  5 


local  orthogonal 
expression 


Option  No.  6 


low-resolution 
amplitude  set 


Option  No.  7 


Option  Nos.  1, 3,  and  4  produce  a  set  of  A(m,j).  Option  Nos.  5.  6.  and  7 
use  an  amplitude  matrix  as  input  and  produce  time  series  or 
algorithms  to  lead  to  time  series.  Option  2  is  a  special  preprocessor 
which  must  precede  the  execution  of  Option  Nos.  3  or  4. 


Figure  2.  SIMBAT  module  relationship . 
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Any  of  the  basic  inputs  (1),  (2,3),  or  (2,4)  can  be  combined  with 
any  of.  the  basic  outputs  (5),  (6),  or  (7).  Typical  example  runs  are 
shown  in  Appendix  A  and  as  files  on  the  accompanying  diskettes.  The 
fundamental  output  for  Option  No.  6  is  given  on  the  diskette  as  file 
MAST6.DAT.  Option  No.  6  requires  further  theoretical  discussion. 
Similarly,  the  concepts  of  conditional  probability  deserve  an  expanded 
exposition. 

7.0  ORTHOGONAL  EXPANSION  (OPTION  6) 

The  shear  mass  of  computations  required  to  compute  forces  at  many 
load  points  in  a  moving  structure  subject  to  wave  actions  is  a  major 
problem  in  operating  with  either  conditional  or  unconditional 
simulations  of  wave  kinematics.  One  approach  is  to  try  to  reduce  the 
number  of  waves  required  to  produce  (approximately)  the  same  wave  train. 
This  is  provided  by  Option  No.  7.  Another  approach  is  to  summarize  the 
local  variations  of  the  wave  kinematics  in  some  sort  of  additive 
function  system.  Option  No.  6  is  based  on  Legendre  and  shifted  Legendre 


orthogonal  polynomials.  Let 

pn(x) 

i  i  i  n 

*  a  +  a.x  +  . . .  +  a  x 
o  1  n 

(28) 

p*(x) 

*  a*  +  a*x  +  . . .  +  a*xn 
o  1  n 

(29) 

be  the  Legendre  and  shifted  Legendre  othogonal  polynomials  of  order  n. 
The  coefficients  are  selected  so  that 

1 

1  p  ( x)  p  ( x)  dx  = 

)  m  n 

-1 

1 

/  pI(x)  pn(x)  dx  = 

0 


0  ,  if  m  M 

2 


2n+l 


,  if  m  =  n 


2n+l 


,  if  m  f  n 
,  if  m  =  n 


(30) 


(31) 
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The  first  several  Legendre  polynomials  are 
Pq(x)  =  1 
Px(x)  =  X 
p2(x)  =  (3x2-1)/2 
p3(x)  =  (5x3-3x)/2 
p4(x)  =  (35x4-30x2+3)/8 
p5(x)  =  (63x5-70x3+15x)/8 

The  SIMBAT  programs  incorporate  Legendre  polynomials  up  to  order  12. 

Most  approximations  to  wave  kinematics  will  not  require  orders  greater 
than  5. 

The  use  of  orthogonal  polynomials  to  approximate  an  arbitrary 
function,  g(x),  defined  on  (-1,1)  can  be  illustrated  as  follows. 

Suppose  the  approximation  to  be  used  is 

N 

g(x)  Z  \  an  Pn(x)  (33) 

n=0 

The  coefficients  afl  are  chosen  from  a  "least-squares"  criterion.  Let 

a  be  those  values  which  minimize 
n 

i  r  n 

Q  =  /  *(X)  •  l  an  Pn(x) 

- 1  [  n=0 

Then 

i  r  n 

I  2  g(x)  -  J  an  Pn(x)  Pk(x)  dx  (35) 

1  n=0 

1  N  1 

i  (a?)  “  ’  /  g(x)  pk(x)  dx  +  J  an  /  Pn(x)  Pk(x)  dx 

*  -1  n=0  - 1 


1S_  =  . 

3ak 


(36) 


Q  is  at  an  extreme  if  3Q/3a^  =  0  for  all  k=0, 1,2, . . . ,N. 

N  1  1 

I  *n  j  pn^x)  pk(x^  dx  =  /  pk(x)  dx 

n=0  -1  -1 

But  by  the  orthogonality  relation,  this  further  reduces 

1 

\(sn)  *  /  8(x)  pk(xl  dx 

-1 


This  reduces  to 

(37) 

to 

(38) 


1 

\  =  (^P)  /  *<*)  Pk(x)  dx 

-1 

The  similar  development  for  shifted  Legendre  polynomials  gives 

1 

a*  =  ( 2k+l )  j  g(x)  p*(x)  dx 

0 


(39) 


(40) 


How  can  these  relations  be  applied  to  ocean  wave  kinematics?  The  essen¬ 
tial  cannonical  form  for  a  linear  wave  property  is  given  in  Equation  24. 
It  should  be  noted  that  in  every  case,  G  (z)  is  either  1.0  (water  level 

IT! 

elevation)  or  is  of  the  form 


-k  z  -k  ( 2d-z) 

in  ,  m 

e  +  Sj  e 


1  +  s2  e 


-2k  d 
m 


(41) 


where  k  is  the  wave  number.  Thus,  the  general  wave  property,  p(n  At) 

ID 

can  be  expressed  from  Equation  10  as 


N-l 


p(n  At)  =  Y  c„ 


i2irmn/N 


(42) 


m=0 


where  C  is  the  FFT  coefficient  given  by 

(Q 


j 

l 


A  ,  e 
mj 


•i&  k  (x  cos(9-0  )  +  y  cos(9-9  )} 
o  m  x  1  7 


(43) 
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for  sea  surface  elevations,  and 


C 

m 


i  { 

)  A_ ,  T_  H,  L 


‘kmz  +  -km(2d-z)l 
s  +  Sj  e 


>  ,  4.  it  , 

mj  m  j 


(  -k  d) 

{*  *  S2  *  I 


j  =  l 


•iB  k  (x  cos(8-0  )  +  y  cos(8-8  )) 
o  in  x  y 


(44) 


for  the  other  wave  properties.  These  can  be  expressed  as  a  sum  of 
products  of  separate  functions  of  x,  y,  and  z  as 


Sea  Surface 


C  =  \  A  ,  e 

m  £  raj 

J=1 

Other  Wave  Properties 

£  A  .  T  H.  f 

*  I  -JU— "~ih  { 

J»1  1  +  a.  a  1,1 


iB  k  x  cos(8-0  )  iB  k  y  cos(8-0  ) 
o  m  x  o  m  y 

e 


-k  z  -k  (2d-z) 

in  .  m 

e  +3^ 


(45) 


-iB  k  x  cos(0-8  )  -iB  k  y  cos(0-8  ) 
o  m  x  o  nr  y 

•  e  e  J 

Suppose  it  is  desired  to  obtain  a  good  representation  locally  in  the 
vicinity  of  the  structure.  Consider  the  volume 

x  -D,  £  x  £  x  +D, 
o  1  o  1 


(46) 


Vd2  *  y  *  yo+D2 


(47) 


0  i  z  i  d 


(Here  z  has  been  taken  as  positive  downwards  and  zero  at  mean  water 
level . ) 
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It  is  natural  to  scale  the  function  as 


Sea  Surface 


:  =  j  a  . 

m  L  &  j 

j=l 


,x-x  > 

■iB  k  x  cos( 8-8  )  -iB  k  D,  cos(9-9 

orao  x  o  m  1\  /  > 


(yy  qV 

-iB  k  y  cos( 9-0  )  -ifl  k  D„  -r-2  cos(9-0  ) 

o  nro  y  o  m  2\  D_  /  y 


Other  Wave  Properties 


r  A  .  T  H,  f  -k  z  -k  (2d-z)! 

:  »  \  _ "J.  .  ■  J  ie  m  +  s  e  m  1 

Z  -2k  d  ie  S1  6  i 


j=l  1  +  s2  e 


-iB  k  x  cos ( 0-9  )  -iB  k  D. (-pj-2)  cos(0-8  ) 
o  m  o  x  o  m  1\  D.  I  x 

e  e  1 


y  -  y 

•iMvn  cos(8-8  )  -iB  k  cos( 9-0  ) 

o  m  o  y  o  m  2\  D n  f  y 


x  «x 

-iB  k  D  2)  cos( 8-8  )  N 

°  ®  D,  /  x  r  ,x-x  . 

-  I  •.  4af) 

_ A  ^ 


-iB  k  D.f-r-2)  cos(0-0  )  N 

° » 2\  °2  /  y  .  x  a  ryo\ 
5  -  Z  aB  pb['dt) 


-k  z 

=  2  c*  pr(e  °  > 
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where  kQ  is  a  selected  single  reference  wave  number.  For  many  applica¬ 
tions,  the  second  term 

-k  (2d-*) 
m 

6 

is  negligible  because  depth  is  large.  For  the  moment  suppose  that  this 
second  term  can  be  ignored.  It  will  be  reintroduced  later.  Then  if 

u  =  (x-xq)/D1 

v  =  (y-yQ)/D2  (53) 

-k  2 
o 


Sea  Surface 


C  =  \  A  .  e 

m  L 

J-l 


-10  k  [x  cos(8-0  )  +  y  cos(8-0  )] 
a  ml  o  x  Jo  y  ‘ 


N  N 

l  \  P«(u)  •  l  b0  p0(v) 
0=0 


o=0 


(54) 


Other  Wave  Properties 
J 


A  ,  T_  Hjjs 


=  J  nm]  m 
j=l 


-i0  k  (x  cos(8-8  )  +  y  cos(8-8  )] 
o  ml  o  x  o  y  ' 


1  +  s2  e 


-2k  d 
m 


N 


N 


N 


-k  z 
o 


l  aa  pa(u)  '  l  b0  p0(v)  -  \  CT  pI(e  ~  ) 


o=o  0=0  r=o 

If  this  is  substituted  into  Equation  42 
N  N  N 


p(n  At)  =  J  J  J  Ba  0  T^n  Pa^U^  p0^v^  pJ^e  °  ^  (55) 


o=0  0=0  T=0 


(54) 
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where  (sea  surface  case) 


B„  ft  r(n  At) 
S,D,i 


N-l 


I 

1=1 


a 

a 


b„  A  . 
3  mj 


-i0  k  x  cos(0-0  ) 
o  m  o  x 

e 


-i0  k  y  cos(0-0  ) 
o  nro  y 


ei2irmn/N 


(56) 


This  is  an  FFT  of  the  quantity  within  {  j  for  each  a,  0,  and  T  combine 
tion. 

The  other  wave  properties  have  a  similar  expression  with 


N-l 


!a,e,r<n  st)  ■  l 


m=0 


J 

l 

j  =  l 


a  b.  cr  A  .  T  H. 
a  p  I  mj  m  j 


-iB  k  [x  cos(0-8  )  +  y  cos(0-0  )] 
o  ml  o  x  Jo  y  ' 


1  +  s2  e 


-2k  d 
m 


i2irmn/N 


.  (57) 

At  a  given  time,  n  At,  many  of  the  B  R  ,  for  a  given  wave  property, 
are  negligible.  After  all  the  region  and  yo±D2  is  relatively 

small  relative  to  the  wave  lengths.  Hence  low  crder  polynomials  are  all 
that  are  required  in  order  to  represent  the  variation  over  the  horizontal 
region.  The  variation  vertically  is  more  or  less  exponentially  atten¬ 
uated  with  depth,  so  a  polynomial  in 

-k  z 
o 
e 

should  only  need  relatively  low  order. 

Hence,  at  a  given  time,  only  a  few  B  -  v  will  be  needed  to  repre- 

®  i  • 

sent  the  wave  property.  The  particular  coefficient  needed  may,  however, 
be  different  from  one  time  step  to  another. 
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Up  to  here  the  actual  computation  of  a^,  b^,  and  c^.  has  been  not 
explicitly  stated.  From  the  definition  of  orthogonal  polynomials 


a 

a 


b 


B 


2g+l 

2 


2B+1 

2 


-ik  0  D,  cos(0-0  )u 

e  p  (u) 

<x 


-i 


l 

-i 


-ik  B  D-  cos(0-0  )v 

m  o  ^  v  /  \  j 

e  7  p^Cv)  dv 


,  k  /k 

=  (2J+1)  I  Wj"  °  Pj(wx)  dw 


(58) 


(59) 


(60) 


8.0  IMPLEMENTATION  IN  CODE 

The  a  ,  bB ,  and  cv,  which  are  functions  of  m  and  j,  are  computed 
a  p  J 

and  combined  as  given  in  Equation  56  and  57  and  then  Fourier  transform 
with  the  fast  Fourier  transform  to  develop  B  ,  _(n  it)  for  each  wave 
property.  These  are  sorted  in  order  of  absolute  value  at  each  time  step 
and  listed  in  a  sequential  file. 

The  coefficients  are  listed  in  the  file  in  integer  form  with  the 
last  3  digits  giving  an  index  which  may  be  used  to  determine  the  orders 
o,  B,  and  1  for  that  coefficient.  Thus,  the  coefficient  with  value 
xxx. xxx  is  listed  as  the  integer  xxxxxxyyy  where  yyy  is  the  order 
designator.  A  matrix 

LSTXYZ  (yyy, 1)  =  a 

LSTXYZ  (yyy, 2)  =  B 

LSTXYZ  (yyy, 3)  =  J 

gives  the  order  associated  with  each  yyy  value.  The  integers  are  ranked 
in  order  of  decreasing  absolute  value.  Thus  the  user  can  compute  with 
the  much  abbreviated  list  of  coefficients  needed  at  that  time  step  to 
represent  the  wave  property  within  the  local  region. 
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9.0  OTHER  DEPTH  TERM 


Let 

-k  (2d -z) 
o 

w2  =  e 

Then  an  exactly  similar  expansion  with  the  same  coefficients  can  be 
developed.  The  resulting  representation  of  the  wave  property  is 

N  N  N 

p(n  At)  =  }  l  l  V 3,y(n  At)  Pa(u)  Pb(v) 

a=o  b=o  r=o 


’  -k  (2d-z)  , 

)  +  Sl  p*(e  ) 


Note:  The  Module  SIM6  is  still  under  testing  and  may  be  changed  further 
as  the  study  continues. 


Appendix  A 
EXAMPLE  RUNS 


OPTION  NO.  1  CONSOLE  LISTING 
TO  BE  FOLLOWED  BY  OPTION  NO.  5 


(■**** 

PROGRAM  SIMULATES  WAVE  PROPERTIES  BY  FREQUENCY  DOMAIN  METHOD 
FOR  EITHER  UNCONDITIONAL  OR  CONDITIONAL  SIMULATIONS. 
(WRITTEN  BY  LEON  BORGMAN,  LARAMIE.  WYOMING) 


PLEASE  KEY  RETURN  TO  CONTINUE 


DO  YOU  WISH  TO  READ  OUTPUT  FROM  PREVIOUS  RUNS 
OF  OPTIONS  41  OR  42?  (Y  *  YES.  N  =  NO 
N 

#******-«•**■**■*  *+****♦♦*  **************  ************* 

ENTER  A  50-CHARACTER  TITLE  FOR  THE  TIME  SERIES  DATA. 


5/38/88  3:07  AM  CONSOLE  LIST  EXAMPLE 

♦h*************-#-#*-*-*-***-*-*  ■*■*•*■*■**■****•**■#*•****■*■**■»■■** 
ENTER  A  Si)— CHARACTER  SUMMARY  OF  THE  TIME  SERIES 
DOCUMENTATION  FOR  FUTURE  REFERENCE. 


THE  OUTPUT  IS  STORED  ON  L: LIST1.DAT 
OPTIONS  ACTIVE  IN  THIS  SUBPROGRAM: 

0.  EXIT  PROGRAM 

1.  PRODUCE  COMPLEX  AMPLITUDES  FOR  UNCONDITIONAL  SIMULATIONS 

2.  PRE-PROCESS  DATA  TO  PRODUCE  NECESSARY  INPUT  FOR  A 

CONDITIONAL  SIMULATION 

3.  SIMULATE  TIME  SERIES  FROM  AMPLITUDES  (USES  OUTPUT  FROM 

STEP  #1.  #3.  OR  #4.  ) 

8.  HELP 


(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 

i  OPTION  NO. 

ENTER  SPECTRAL— LINE  CUTOFF  FRACTION.  (NOTE:  A  SMAi_u 
VALUE  IS  GOOD.  ONLY  THOSE  LINES  GIVING  VARIANCE  CON "R 1 9.“ I  0 
GREATER  THAN  OR  EQUAL  TO  ( CUTOFF*LA.RGEST  SPECTRAl.  uInE)  ARE 
KE0T 

a.  00001 

*##***♦-***  *******•*•»*>******-***************  ****>** 

ENTER  SEED  INTEGER  FOR  RANDOM  NUMBER  GENERATOR.  1 13 


123436 


ENTER  I OPT: • 

I OPT- l  INDICATES  DSPEC  MATRIX  IS  COMPUTED  AND  THEN 
STORED  FOR  FUTURE  USE  IN  A  USER-SPECIFIED  FILE. 

'  UTi.:  T  \in  T  roT?q  .mqti'y  »c  o^ro n  —  m  a 


**#****♦♦**♦****•*****■**■**■**•**•«■■*•*  ******************** 
1 

ENTER  FILE  NAME  TO  WHICH  THE  DATA  IS  TO  BE  STORED. 
JsDSPECl.DAT 

#•*■»■»•***-*■*♦-***■»•********■****■)*•***********#***•»***** 

ENTER  A  30-CHARACTSR  TITLE  FOR  DRSPEC.  MATRIX 


EXAMPLE  DIR.  SPECTRA  NO.  1 

MODE  NUMBER  AT  THIS  STEP  =  1 

DO  YOU  WANT  TO  ENTER  ANOTHER  MODE? 

ENTER  N-NO  IF  NO  MORE  MODES  ARE  “0  BE  ENTERED. 
OTHERWISE.  ENTER  Y-YES 
WHAT  IS  YOUR  CHOICE? 

Y 

ENTER  LAMDA  FOR  0-H  SPECTRA.  NOTE:  LAMDA-1.0  FOR 
P-M  SPECTRA  OR  JONSWAP  SPECTRA 
1.  0 

ENTER  MODAL  FREQUENCY  FOR  MODE. 

'll.  1 

ENTER  TOTAL  VARIANCE  OF  MODE. 

VARIANCE  DIMENSIONS  DETERMINE  SPECTRA  DIMENSIONS 
DR. SPECT. DENSITY  WIU_  BE  IN  UNITS  OF 
LENGTH**^  PER  (HERTZ-RADIAN) 

1012.  0 


PRINC. DIR. -PRINC.  DIR.  CONST.  + 

PR INC. DIR.  SLOPE* (FREQ.  -MODAL  FREQ.  ) 
DIMENSIONS:  PRINC. DIR. CONST.  IN  NAV.  DEGREES 

princ.  dir. slope  in  degrees— sec. 

ENTER  PRINCIPAL  DIRECTION  CONSTANT  FOR  MODE. 

0.  0 

ENTER  PRINCIPAL.  DIRECTION  SLOPE  FOR  MODE. 

0.  0 


3PRD. STD. DEV.  -SPRD.  STD.  DEV.  CONST.  + 

SPRD. STD. DEV.  SLOPE* (FREQ. -mODAu  FR£Q.  ) 

DIMENSIONS:  SPRD. STD. DEV. CONST.  IN  NAV.  DEGREES 
SPRD. STD. DEV. SLOPE  IN  DEGREES-SEC. 

ENTER  SPREADING  STD.  DEV.  CONSTANT  FOR  MODE. 

30.  0 

ENTER  SPREADING  STD.  DEV.  SLOPE  FOR  MODE. 

0.0 

SPECTRA  PARAMETERS:  LAMDA.  F0.  VAR-  1.00000  .10000  100. 

SPREAD  PARAMETERS:  THET0. THET1, SIG0. SIG1-  .00000  .00000 

ARE  THESE  THE  VALUES  YOU  WANTED?  IF  NOT,  ENTER  N-NO 
AND  RE-ENTER  PARAMETERS  FOR  THIS  MODE. 

OTHERWISE  ENTER  Y-YES 
WHAT  IS  YOUR  CHOICE? 

Y 


MODE  NUMBER  AT  THIS  STEP  =  £ 

DO  YOU  WANT  TO  ENTER  ANOTHER  MODE? 

ENTER  N-NO  IF  NO  MORE  MODES  ARE  TO  BE  ENTERED. 
OTHERWISE.  ENTER  Y-YES 
WHAT  IS  YOUR  CHOICE? 

N 


AMPLITUDE  RANDOMNESS  MENU 

1.  RANDOM  PHASE.  BUT  WAVE  AMPLITUDE  DETERMINISTIC  AND 
CONSTRAINED  TO  BE  EQUAL  TO  E.  0*SQRT (S (F.  THETA) ) 
E.  RANDOM  PHASE  AND  RANDOM  AMPLITUDE 

please  enter  your  choice 


POINT  B  REACHED-' 

POINT  C  REACHED 
POINT  D  REACHED 

NUMBER  OF  DEGREES  OF  FREEDOM  =  1737 
OPTIONS  ACTIVE  IN  THIS  SUBPROGRAM: 

3.  EXIT  PROGRAM 

1.  PRODUCE  COMPLEX  AMPLITUDES  FOR  UNCONDITIONAL  SIMULA" IONS 
£.  PRE-PROCESS  DATA  TO  PRODUCE  NECESSARY  INPUT  FOR  A 
CONDITIONAL  SIMULATION 

3.  SIMULATE  TIME  SERIES  FROM  AMPLITUDES  -'USES  OUTPUT  FROM 
STEP  #1.  #3.  OR  #4.) 
a.  HELP 


•»■***** 


(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 


3 

DO  YOU  WANT  TO  SAVE  THE  OUTPUT  FROM  LAST  OPTION 
STEP  ON  A  SEPARATE  USER  SPECIFIED  FILE  ?  Y  OR  N 
Y 

PLEASE  ENTER  THE  FILE  NAME  FOR  STORAGE 

J :  CUT  1 .  DAT  OUTPUT  FILE 

Programmed  STOP 


The  printer  listing  Is  In  file  UST1.DAT. 

The  output  to  be  used  as  input  to  Option  no.  5  is  stored  in  OUT2.DAT. 


L:SIM12' 


CONSOLE  LIST  FOR  OPTION  NO.  5 
FOLLOWING  OPTION  NO.  1 


PROGRAM  SIMULATES  WAVE  PROPERTIES  5Y  FREQUENCY  DOMAIN  METHOD'S 
FOR  EITHER  UNCONDITIONAL  OR  CONDITIONAL  SIMULATIONS. 

(WRITTEN  BY  LEON  BORGMAN,  LARAMIE,  WYOMING) 


PLEASE  KEY  RETURN  TO  CONTINUE 


DO  YOU  WISH  TO  READ  OUTPUT  FROM  PREVIOUS  RUNS 
OF  OPTIONS  #1  OR  #£?  <Y  =  YES.  N  =  NO 
Y 

WHAT  IS  THE  FILE  NAME  FOR  PREVIOUS  OUTPUT? 

J:OUTl. DAT 

OPTIONS  ACTIVE  IN  THIS  SUBPROGRAM: 

0.  EXIT  PROGRAM 

1.  PRODUCE  COMPLEX  AMPLITUDES  FOR  UNCONDITIONAL  SIMULATIONS 

2.  PRE-PROCESS  DATA  TO  PRODUCE  NECESSARY  INPUT  FOR  A 

CONDITIONAL  SIMULATION 

5.  SIMULATE  TIME  SERIES  FROM  AMPLITUDES  (USES  OUTPUT  "ROM 
STEP  #1,  #3.  OR  # 4.  ) 

8.  HELP 

.ft**.****  ***##**•******#***##*#*##•*#•****#■*•**•*••**•*••►■*■*•** 

(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 

S  OPTION  NO. 

ENTER  NOISE  TO  SIGNAL  RATIO.  AS  RATIO  OF  NOISE 

STANDARD  DEVIATION  TO  SIGNAL  STANDARD  DEVIATION 

0.  0 

CHANNEL  NUMBER:  1 

PLEASE  ENTER  MEAN  VALUE  DESIRED  FOR  THIS  CHANNEL 

0.  0 

CHANNEL  NUMBER:  £ 

PLEASE  ENTER  MEAN  VALUE  DESIRED  FOR  THIS  ChANNE^. 

0.  0 

CHANNEL  NUMBER:  3 

PLEASE  ENTER  MEAN  VALLE  DESIRED  FOR  THIS  CHANNEL 

0.  0 

ENTER  TIME  STEP  FOR  START  OF-  LIST  AND  RETURN 
1 

ENTER  TIME  STEP  AT  TERMINATION  OF  LIST  AND  RETURN 
£00 

ENTER  0  IF  GRAPHICAL  OUTPUT  IS  DESIRED 
KEY  1  IF  NUMERICAL  LIST  OF  TIME  SERIES  IS  WANTED 
0 

NUMBER  OF  DEGREES  OF  FREEDOM  =  0 

'ISTfnva  Of~— T(,.‘C  TV  TL.T8  «t  1“  OOnfSQQTVt  • 


£.  PRE-PROCESS  DATA  TO  PRODUCE  NECESSPRY  INPUT  FOR  P 
CONDI TIONPL  SIMULATION 

5.  SIMULATE  TIME  SERIES  FROM  AMPLITUDES  (USES  OUTPUT  FROM 
STEP  #1.  #3.  OR  #4.  ) 

3.  HELP 

************************************************** 

(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 


ENTER  NOISE  TO  SIGNAL  RATIO.  AS  RATIO  OF  NOISE 

STANDARD  DEVIATION  TO  SIGNAL  STANDARD  DEVIATION 

0.  0 

CHANNEL  NUMBER: 

PLEASE  ENTER  MEAN  VALUE  DESIRED  FOR  THIS  CHANNEL 

0.  0 

CHANNEL  NUMBER:  £ 

PLEASE  ENTER  MEAN  VALUE  DESIRED  FOR  THIS  CHANNEL 

0.  0 

CHANNEL  NUMBER:  3 

PLEASE  ENTER  MEAN  VALLE  DESIRED  FOR  THIS  CHANNEL 

0.  0 

ENTER  TIME  STEP  FOR  START  OF  i_IST  AND  RETURN 
1 

ENTER  TIME  STEP  AT  TERMINATION  OF  LIST  AND  RETURN 
£00 

ENTER  0  IF  GRAPHICAL  OUTPUT  IS  DESIRED 
KEY  I  IF  NUMERICAL.  LIST  OF  “I ME  SERIES  IS  WANTED 
1 

NUMBER  OF  DEGREES  OF  '"REEDOM  ^  0 
OPTIONS  ACTIVE  IN  “HIS  SUBPROGRAM: 

0.  EXIT  PROGRAM 

1.  PRODUCE  COMPLEX  AMPLITUDES  FOR  UNCONDITIONAL  SIMULATION 
S.  PRE-PROCESS  DATA  TO  PRODUCE  NECESSARY  INPUT  FOR  A 
CONDITIONAL  SIMULATION 

5.  SIMULATE  TIME  SERIES  FROM  AMPLITUDES  (USES  OUTPUT  FROM 
STEP  #1.  #3.  OR  #4.  ) 

8.  HELP 

************************************************** 

(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 


DO  YOU  WANT  TO  SAVE  THE  OUTPUT  FROM  LAST  OPTION 
STEP  ON  A  SEPARATE  USER  SPECIFIED  FILE  ?  Y  OR  N 

N  NO  OUTPUT 

Programmed  STOP 


F> 


■»  _:3IM1£5 

CONSOLE  LIST  FOR  PREPROCESSOR 
OPTION  NO.  2 


PROGRAM  SIMULATES  WAVE  PROPERTIES  BY  FREQUENCY  DOMAIN  METHODS 
FOR  EITHER  UNCONDITIONAL  OR  CONDITIONAL  SIMULATIONS. 

(WRITTEN  BY  LEON  BORGMAN.  LARAMIE.  WYOMING) 

PLEASE  KEY  RETURN  TO  CONTINUE 


DO  YOU  WISH  TO  READ  OUTPUT  FROM  PREVIOUS  RUNS 
OF  OPTIONS  *1  OR  #£?  (Y  *  YES.  N  *  NO 
N 

***#*****************•***********.***.*.*****.*.*.****.*. 

ENTER  A  50— CHARACTER  TITLE  FOR  THE  TIME  SERIES  DATA. 


5/33A93  3:£6  PM  CONSOLE  LIST  FOR  OPTION  £ 

*******■**********■*•*■**********■*.**.*.)»..******.**.*.**.*** 

ENTER  A  50-CHARACTER  SUMMARY  OF  THE  TIME  SERIES 
DOCUMENTATION  FOR  FUTURE  REFERENCE. 


-1ST  OUTPUT  STORED  ON  JL  J:LIST£.  DAT 
OPTIONS  ACTIVE  IN  THIS  SUBPROGRAM; 

0.  EXIT  PROGRAM 

1.  PRODUCE  COMPLEX  AMPLITUDES  FOR  UNCONDITIONAL  SIML-ATICNS 
£.  PRE-PROCESS  DATA  TO  PRODUCE  NECESSARY  INPUT  FOR  A 
CONDITIONAL  SIMULATION 

5.  SIMULATE  TIME  SERIES  FROM  AMPLITUDES  (USES  OUTPUT  FROM 
STEP  #1.  #3.  OR  #4.) 

3.  HELP 

#**********#**************************.*.»..*.**.*.*.*.*.**.* 

(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 


£  OPTION  NO. 

ENTER  SPECTRAL-uINE  CUTOFF  FRACTION.  (NOTE:  A  SMAi.__ 

VALUE  IS  GOOD.  ONLY  THOSE  LINES  GIVING  VARIANCE  CJnTR I  BUT I -N 
GREATER  THAN  OR  EQUAL  TO  ( CUTGFF*LARGES“  SPECTRAL  LINF)  A  •?" 
KEPT 


0.  000001 


•»****#** 


ENTER  SEED  INTEGER  FOR  RANDOM  NUMBER  GENERATOR.  110 


1337334 
ENTER  I OPT; 

I  OPT*  1  INDICATES  DSPEC  MATRIX  IS  COMPUTED  AND  T.-.EN 
STORED  FOR  FUTURE  USE  IN  A  USER-SPECIFIED  FILE. 


i>V'- 


*****************  ******  *****  ***  *******  ************** 

5 


POINT  V  REACHED 

ENTER  =  ILE  NAME  FROM  WHERE  THE  DATA  IS  TO  BE  READ. 
JsDSPECl. DAT 
POINT  W  REACHED 
POINT  X  REACHED 
POINT  Y  REACHED 
POINT  Z  REACHED 

ft***********'*******#*****'*"*************#****'***** 

AMPLITUDE  RANDOMNESS  MENU 

1.  RANDOM  PHASE,  BUT  WAVE  AMPLITUDE  DETERMINISTIC  AND 
CONSTRAINED  -rC  BE  EOUAu  TO  2. 0*SQRT <S CF.  THETA) ) 
£.  RANDOM  PHASE  AND  RANDOM  AMPLITUDE 
PLEASE  ENTER  YOUR  CHOICE 


POINT  B  REACHED 
POINT  C  REACHED 
POINT  D  REACHED 

NUMBER  OF  DECREES  OF  “REEDOM  =  1729 
OPTIONS  ACTIVE  IN  THIS  SUBPROGRAM: 

3.  EXIT  PROGRAM 

1.  PRODUCE  COMPLEX  AMPLITUDES  FOR  UNCONDITIONAL  SIMULATIONS 

2.  PRE-PROCESS  DATA  TO  PRODUCE  NECESSARY  INPUT  FOR  A 

CONDITIONAL  SIMULATION 

S.  SIMULATE  TIME  SERIES  FROM  AMPLITUDES  (LSES  OUTPUT  FROM 
STEP  #1.  #3.  OR  #4.  ) 

9.  HELP 

**  *******************  ***************************** 

{PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 


DO  YOU  WANT  TO  SAVE  THE  OUTPUT  FROM  LAST  OPTION 
3-EP  ON  A  SEPARATE  USER  SPECIFIED  FILE  ?  Y  OR  N 
J : CUT 2. DAT 

Programmed  STOP 

r  > 


Store  output  on  OUT2.DAT. 
Printer  output  is  shown  on  UST2.DAT. 


CONSOLE  UST 
OPTION  NO.  3 


*** •*■>*■!<•**  »•#■«•*■)»••»*  *  *-**  ***Hir*-«r  )»■****•«■**■*•  ►*  **■»  4 •*■**•*•*  *•* 

PROGRAM  SZ»siiJwfiTES  WAVE  — 1 'RGPEATIES  U Y  FRE'JlENwV  Ol;MA - N  ,y El  .  — 0 

for  either  unconditional  or  conditional  situ-atio-ms. 

(WRITTEN  BY  _SQN  BGR3MAN.  LARAMIE,  WYOMING) 

PLEASE  KEY  RETURN  TO  CONTINUE 

***********  ****************  ***■**■*■**•*■***•*.*■*  •*■*■*•*■#■***■*•*■*>•**  •*•».• 


DO  YOU  rtlSn  TO  READ  OUTPUT  .-ROM  PREVIOUS  R_NS 
0“' OPTIONS  OR  ^2?  (Y  =  YES.  N  =  NC 
Y 

/•-•AT  IS  "  -E  rI_E  nA-^E  .’-OR  -  REV Cl5  ClTPcT? 

J  :  CUTE.  DAT  INPUT 
TITLE:  5/23/03 
DOCUMENTATION:  6: 20  ="« 

ISEED-  12345677  ILHUZ*  2  n-  502 

OPTIONS  ACTIVE  IN  7. R IS  SUBPROSHAM: 

0.  EXIT  PROGRAM 

3.  DEVELOP  COMPLEX-VALUED  AMPLITUDES.  CONDITIONED  ON  GIVEN 
INPUT.  WHERE  COND I ~ I GN I NG  INPUT  IS  _ESS  THAN  N. 

(USES  OUTPUT  "ROM  3  rEP  L2.  ) 

8.  HELP 

**  **■»■* *****#*♦********#*******♦»**•»*+***♦♦+*♦**+«.  * 

(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 


3  OPTION  NO. 

MENU  FOR  OPTION  ±>3  CONDITIONING  I'PUY; 

1.  READ  THE  CONDITIONING  TIME  SERIES  INT£.RVA_ 

2.  INPUT  FROM  CONSOLE  A  SINGLE  POIN"  C-vECK  COM 


1 

THE  GIVEN  TIME  SERIES  IS  TO  BE  READ  FROM 


J:GTS3.  DAT 

AT  A.  XC(1.  1,  i>*+l. 387074E-H32 
AT  B,  XC ( 1 .  1 .  1 >  ~+ 1 . 387 ?74£f 02 


ITER-1 

ITER-2 

ITER-3 

ITER-4 

ITER-5 

ITER-6 

ITER-7 

ITER-Q 

I TER *9 

ITER-10 

ITER- 11 

"COs'  i 


ERROR— 1-2.  87978591 00633 1  E-001 
ERROR— *-1.  49614657735246E-001 
ERROR— t-S.  3367809756308 IE-002 
ERR0R-+6. 2565Q778317937E-002 
ERROR— *-4.  79428 186085693E-002 
ERROR— *-3.  986557227 17029E-002 
ERROR— 1-4.  07936739941 164E-002 
ERROR— *-4.  8252003641 025 IE-002 
ERROR— f 4. 0 1 262099626422E— 002 
ERROR— +3. 7075467994033 IE-002 
ERROR— *-3.  61 28806 1379002E-002 

TOO  no**-?  70,7|07A7?4  <  Q  777C_i7|.7);r. 


-'ROM 


rr. 


ITER=15 
I  TER® 16 
I TER=1 7 
ITER® 13 
I-ER®13 
*7ER®20 

I“ER=22 


ERRGR=-4.  315727206652996-002 
6RR0R=+3. 6527660 73345 11 E-002 
ERROR®-*. 05362967637325E-002 
ERRQR®-4.  3041 1327745454E-002 
ER.R0R=+4. 2302564824 1650E-002 
ZRRCR®+3. 22563350 135396- -002 
ZRRC3R®+4.  523974605323202-002 
ERRCR=+5. 3067S454774573E-002 
ERR0R=-6.  322746 7621 0544E- 002 
2,R,RGR=+8. 509675 19543353E-002 
ERROR®- 1. 311356241472206-001 
2RR0R®+1. 776503371 12263E-001 
ERRCR®*3.  4 1 351 364S35S63E-002 
ER!RGR=-4. 44513361 145773E-002 
ERRCR«-7. 1 5560563985SS3E— 002 
ERROR®-*. 29335306096 395E-002 
IR.RCR=+4. 51 4495349223 17E-002 
ERROR®-*. 1 9 1 3239 161 9607 2— 002 
ERRCR=+3. 42003073337953E-002 
ERRO 7**4. 777794130333736-002 
=  R  ■*’ 3 R®  — 2.  3ooob365£ 32  o  4*2  o— 008 
”  R  ROR®— *.  1  o75304fei33£  1  332 — *002 
ERROR®— 5.  *  '-372632905 1 4*£— 002 
2  R  RG :R*  -5 .  464*4  ;  032220246-  002 
oR^GR®— 6*  33353233327641=.— 002 
=  RRG  234054o  1  S'! 40 12o—002 
I’RCR-'a. 067366397222036-002 
E RRC R®-.-*. 39337506333 04SE-002 
I  R  RC  R— **  .  06 1 3 1 9  =4=*  3 038 z, —  7 0 = 


I  R  R,:  !-  .  02463503219*331  =—001 

ERROR® -3. 905330356644306-002 
6  REG  R— 1 .  1  *5259022263376-00 1 
6RRCR—2.  043052720434036-002 
ERRGR=-2. 60553536376 ; iSE-002 
o  R  .R--  E  —  *■  t  *  1 14 o*  /si  1371 3bo2-002 
IRRE  R®-h. 3936 130330*4 136-002 
1".  RCR®  -2.  2:56376663337 366-302 
■-  .R _  R  — — 8 .  1 3  w  ‘a 5 9 4 3= 5 o  7  o  o  2 — 0  0  2 
-:ER.:E®~3.  373773337895236-002 
E  RRGR®  *-3.  6873 1 983806  3726— 002 
ERGR»— 8.  358o39o6056o252-008 
6 E.R0R®-2. 342238732122436-002 
6  ERCiE®- 1 -  54923490 1 327496-002 
ERROR®-!. 653733954793336-002 
ERROR®— 5.  662:3367725365 1 6-  302 
6REGR=-3. 63765590 127990E-002 
IRRC'-R®  -i.  57571394050733 


G  E  E«  -2.  1 32  7695  ill  464 1 « 
..  R.RGR-4.6.  S  10025.73749505 
6  RRC  E®-*.  34095337*31 033: 


1  E  ;c  -4.  07719354531 7336-00 
6  R.RGR®  i-3.  667 1 07 9*6,044066-00 
E  R  EGR®+3.  53025320207473E-00 
ERROR®**. 932297934445236-00 
ERROR®-*. 297 17308033476E-00 
E.R.RC  E®  -3.  0363337303531 7  E-00 
TRRCR---+2.  14243542 172331  E-00 
ERROR®-:. 094708 13600263E-00 
ERROR®-: .  0:0452704605772-00 
ERROR®- 1 . 60*9.7317293156-00 


i  nj  fij  OJ  ,>j  nj  fij  ftj  i;.j  OJ  OJ  OJ  OJ  > 


m  H  il  .'J  Al  .-J 


-/■j  V  -  _  .j  . -■  v  _  ,  ~r.  — w  -  -  -  -  \ji7~  — --Vj  Lj  — '  i 

S~E“‘  L\  A  GEP3RA”£  USES  3=>2C:-IZD  -ILE  ?  V  QR  V 

i 

*'„£.-SE  E\— E3  "  -E  ”-_Z  3”CRPGE 

- !  u-“3.  -AT  OUTPUT  FILE 

-rorra.-ornea  £TCP 

r-> 


til  tij  (51I11  15  r j  11  15  1.1  i'i  >5  111  -■»  Ml  r.|  lit 


CONSOLE  UST  FOR  OPTION  NO.  5 
WITH 

OUT3.DAT  AS  INPUT 


*  *.4  **■*  ****<•**•*■**#>*•**■*  **■•*■*•* *-*  *-***-*^  *•**■►*■**  *♦*-*■*-*•*  ♦ 


-a ;r  3 1  Y  JL.fiT<33  vs^Vili  v RGP!i*  r3T 1 33  3Y  r~3EQiJ£NCY  D G !yl PIN 
t-R  £;--.£*•  wMCONDITID.NAU  OR  C DND I T I GNflL  SIMUuATCnS. 

;  NR  I  **70  -  BY  LSCN  BCRSi*tP«v,  LARAMIE.  WYOMING) 


PLEASE  KEY  RETURN  70  CONTINUE 

*  *  *■*#*■+*#+  +  +  **•* •*•*•*•*•**•»' •**  *«*■#•  * -*■•*•***■«■  *•**■*■  •*****■*■  ■*•*■*••*  •**•*•**•■*•#••»•■*•  «• 


/  —  t  —  _  *v  .  o  — . 


“0  ™'ERC  OJT-'j”  ”7  Cm 

-  OR  *:cT'  •.  Y  VE3,  N 


-•73YI0-3  RUNS 
_  , 


■t  r> 


**  *'*■1  •  *  _  -0.  -  ^  _  -  ■*.  -  **  —  L.*  L 

-VE-T.  3  47  INPUT 

0 ~  I  0 '•  3  0  .  Y  E  I '■  '*'  •10  3  _Z‘  ”  '•C37A  m : 

a.  ox: 7  program 

.*..  -'RCDuDE  CCMP.Ea  AmPl.”UDES  -OR  UNCONDITIONAL  SiMLu 
0.  PRi— PS0CSS3  DATA  — 0  PRODUCE  NECESSARY  INA'wT  ”07  A 
OO.'iD*  ■■  iw-inL  S..TJ-A  luN 

5.  SIMULA”  “IMS  SERIES  r ROY  AMPLITUDES  <lS£S  OUTPUT 

3  i  if-*  *1.  Tt  o  ■  — :  R  i*A.) 


-  RUM 


•♦**■*■«  *-*►**  *-****■*  +  ***>**■**•*••****■*♦*  •»***■*•***•*•**•*■*#■*#* 

•: -TEASE  ENTER  YOUR  CHOIOE  AND  KEY  RETURN) 

OPTION  NO. 

NT  SR  NOISE  TO  S:S\Aw  RATO,  AS  RATIO  OR  NOISE 

STANDARD  DEVIATION  ”0  SIGNAL  STANDARD  DEVIATION 


L-I* 

<v„vr_ 

NUMBER: 

-t5c 

£ 

enter  mean 

VPLL-Z 

DESI 

RED 

T  -13 

CAA.nnE^ 

I  '*4  1  _ 

'■i U M B- R  : 

.z  -*3«^ 

t'.V’n'I. 

E'T'E.R  mean 

VS-L..Z 

DEE  I 

RED 

r  Z.9 

T  n  1 3 

C,- A  v  n£._ 

Nu  v  BE  R : 

,  3  •*  3  3 

iN  *ER  mEAt-i 

VAi— Lit 

D -5  .1 

RED 

C'”: 

h  1 3 

C.-AnnE- 

.  ,3 

N TER  TIME  STEP  ”CR  S”AR7  0-  lIS”  A\<D  RETURN 

NTR  TIME  3”EP  A”  “ZRMI.NATOIN  0“  LIS”  AND  RETURN 
d& 

v”SR  a  I”  GRAPHICAL  CUT-UT  IS  DESIRED 
KEY  :  I-  NUMERICAL  UIS”  Or  TIME  SERIES  IS 


-'"BE-*  Or  0E3REES  Cr  “REEDC"  *=  3 


WANTED 


5.  ,3..  '‘UL-'i  3.  -  r’—  3  —  “(I— *3  (  LS-3  u«  r"'j  .  ^ 

S“i='  "- ,  *?3.  Cl  J,*4.  ) 


3. 

**»♦*•*•»♦*■•*■*•**»■*  »*-»- *•»*«"»>*.  -*«•*■•«>•■< 

(?_£?3E  ENTER  YCLS  CHC  ICE  «,\D  KEY  1E~ URN ) 


.3 

S3  YC_  *SN'  “C  SAVE  7-E  “lO:*1  _.A3~  G-'~ I cm 


SIM4 


CONSOLE  LIST 
OPTION  NO.  4 


**************************************************** 

PROGRAM  SIMULATES  WAVE  PROPERTIES  BY  FREQUENCY  DOMAIN  METHODS 
FOR  EITHER  UNCONDITIONAL  OR  CONDITIONAL  SIMULATIONS. 

< WRITTEN  EY  LEON  BORGMAN.  LARAMIE.  WYOMING) 

PLEASE  KEY  RETURN  TO  CONTINUE 

****************************  ****************************** 


DO  YOU  WISH  TO  READ  OUTPUT  FROM  PREVIOUS  RUNS 
OF  OPTIONS  #1  OR  #2?  <  Y  -  YES.  N  »  NO 
Y 

WHAT  1 3  T  “  ILE  ’  AmS  ■■  l R  -l:,!iV*iJU'3  O'-i  i  —‘U  t  ? 

JsCL72.DAT  «* -  INPUT 

TITLE:  5/23/S3 
DOCUMENTATION:  6 : 2'2  PM 

ISEED-  1  £345677  ICHUZ-  2  N-  512  NTS* 

OPTIONS  ACTIVE  IN  THIS  SUBPROGRAM : / ’  3.  EXIT  PROGRAM 

4.  DEVELOP  COMPLEX-VALUED  AMPLITUDES.  CONDITIONED  ON  31 V 
INPUT.  WHERE  CONDITIONING  INTERVAL  E0UAL3  N. 

(USES  OUTPUT  FROM  S~EP  42.  ) 

Q.  HELP 

************************************************** 

(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 

4  - - OPTION  NO. 

THE  GIVEN  TIME  SERIES  IS  TO  BE  READ  FROM  WHAT  -:lE? 


J : GTS4. 1 

DAT 

OPTION 

4: 

JM-1 

OPTION 

4: 

JIM-2 

OPTION 

4: 

JIM— 3 

opticn 

4 : 

JIM-4 

OPTION 

4: 

JM— 5 

OPTION 

<* : 

JM-6 

OPTION 

4: 

J  *i=7 

OPTION 

4 : 

jM-a 

OPTION 

4 : 

JIM-S 

OPTION 

4: 

jM=ua 

option 

4: 

JM— 1  1 

OPTION 

4: 

JM-12 

OPTION 

4 : 

JIM- 13 

OPTION 

4: 

JM-  1 4 

OPTION 

4 : 

J  M-  1 5 

OPTION 

4: 

JIM—  1  & 

OPTION 

4 : 

JM—  t  7 

OPTION 

4: 

jM-ia 

OPTION 

4 : 

JM—  1 3 

These  can  easily  be  removed  from 
FORTRAN  code. 

However,  it  le  pleasant  to  see  on 
die  console  that  the  computer 
le  moving  along  In  the  calculations. 


OPTION  4 
OPTION  4 
OPTION  4 
OPTION  4 
OPTION  4 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 


J  1*1=31 
J!*l=3£ 
J*i=33 
JM=34 
JT=35 
JT=36 
JT=37 
JN=38 
J*!=39 
JM=40 
Jy 1=4 1 
JT=4;= 

J'r*  =  44 
J.r=45 
JI*=4S 
JT=47 
Jr  =  48 
jr=49 

Jr  =5i3 


i  i  ?•'  ±  i  **  Z  ■”  *z$ 

OPTION  4:  J*t=89 
OPTION. 4:  JM=90 

OPTION  4:  J’1='3i 
OPTION  4:  JM=9S 
OPTION  4:  J *1=93 
OPTION  JM=94 
OPTION  4:  J  '1=35 
OPTION  4:  JM=36 

OPTION  4:  JM=97 
OPTION  4:  J*=9a 

OPTION  4:  j"M=99 
OPTION  4:  JM=100 
OPTION  4:  JM=101 
OPTION  4s  JM=I0£ 

OPTION  4s  JM=103 
OPTION  4s  J  v  =  1 04 
OPTION  4s  jN=i0S 
OPTION  4s  JT=10£ 

OPTION  4s  JM=> ,1  07 
OPTION  4s  JT=108 
OPTION  4s  09 

OPTION  4s  JM=:10 
OPTION  4s  JM»111 
OPTION  4s  JIH=:i£ 

OPTION  4s  13 

OPTION  4s  .!«=•.  14 
OPTION  4:  J.vi»:  IS 
OPTION  4S  JN=-  1 6 
OPTION  4s  J.>'=i:7 
OPTION  4  s  -rN*r.  I  a 
OPTION  «•:  JM=al9 
OPTION  4s  JN*I£0 
OPTION  4s  JM=I£1 
OPTION  4s  JM»,c;2 
OPTION  4s  JM* 1 23 
OPTION  4s  JN=»I£4 
OPTION  4s  JN=r.£S 
OPTION  4s 
OPTION  4s  JN=t3 7 

OPTIONS  ACTIVE  1%  “-IS  SUBPROGRAM s / 1  0.  EXIT  PROGRA 

4.  DEVELOP  COMPLEX— VALUED  AMPLITUDES,  CCADi~:GNED  0 
INPUT,  WHERE  CONDITIONING  INTERVAL  EOLA-S  N. 
(USES  OUTPUT  FROM  STEP  #£.  ) 

3.  HELP 

a********************************************  ***** 

(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 

0 

DO  YOU  WANT  TO  SAVE  THE  OUTPUT  FROM  LOST  OPTION 
STEP  ON  A  SEPARATE  USER  SPECIFIED  FILE  ?  Y  CR  N 
Y 

please  enter  the  FI_S  NAME  for  S-CRAGE 

J  s  OUT 4.  DAT  OUTPUT  FILE 


Pro a rammed  STOP 


CONSOLE  LIST  FOR  OPTION  NO.  5 
USING  OUT4.DAT  AS  INPUT 


**#********■»■**■»■«•**•«.**•*.**♦******«.*)(.*****#.  ********  *  *  * 

PROGRAM  SIMULATES  WAVE  PROPERTIES  BY  FREQUENCY  DOMAIN  METHOD 
“OR  EITHER  UNCONDITIONAL  OR  CONDITIONAL  SIMULATIONS. 
(WRITTEN  BY  LEON  BORGMAN.  LARAMIE,  WYOMING) 

PLEASE  KEY  RETURN  TO  CONTINUE 

ft**********************  ******-****«- '»**#**•»*****#'***'*-»  .*.***..*.*. 


DO  YOU  WISH  TQ  READ  OUTPUT  FROM  PREVIOUS  RUNS 
OF  OPTIONS  OR  #2?  (Y  =  YES.  N  a  NO 
Y 

WHAT  IS  THE  FILE  NAME  FOR  PREVIOUS  OUTPUT? 

JiOUTA.DAT  INPUT 
OPTIONS  ACTIVE  IN  THIS  SUBPROGRAM: 

0.  EXIT  PROGRAM 

1.  PRODUCE  COMPLEX  AMPLITUDES  FOR  UNCONDITIONAL  SIMULATIONS 

2.  PRE-PROCESS  DATA  TO  PRODUCE  NECESSARY  INPUT  FOR  A 

CONDITIONAL  SIMULATION 

5.  SIMULATE  TIME  SERIES  FROM  AMPLITUDES  (USES  OUTPUT  FROM 
STEP  #1,  #3,  OR  #4. ) 
fl.  HELP 

ft****#*****#*****************.*.*’**  ***«*#  turn- *•*****•»•* 

(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 

5  OPTION  NO. 

ENTER  NOISE  TO  SIGNAL  RATIO.  AS  RATIO  OF  NCISE 

STANDARD  DEVIATION  TO  SIGNAL  STANDARD  DEVIATION 

0.  0 

CHANNEL  NUMBER:  1 

PLEASE  ENTER  MEAN  VALUE  DESIRED  “OR  THIS  CHANNEL 

0.  0 

CHANNEL  NUMBER:  2 

PLEASE  ENTER  MEAN  VALLE  DESIRED  FCR  THTS  CHANNE' 

0.  0 

CHANNEL  NUMBER:  3 

PLEASE  ENTER  MEAN  VALUE  DESIRED  FOR  T,-iIS  CHANNEL 

0.  0 

ENTER  TIME  STEP  FOR  STAR‘r  OF  LIST  AND  RETURN 
1 

ENTER  TIME  STEP  AT  TERMINATION  OF  LIST  AND  RETURN 
200 

ENTER  0  IF  GRAPHICAL  OUTPUT  IS  DESIRED 
KEY  1  IF  NUMERICAL  LIST  OF  TIME  SERIES  IS  WANTED 
1 

NUMBER  OF  DEGREES  OF  “REEDOM  ■  0 


1.  .-‘^wu  uwi  <  w-4wC.-.iji  3_.|_ 

2.  PRE-PROCESS  DA~A  TO  PRODUCE  NECESSARY  INPUT~FOR  A~ . 

CONDITIONAL  SIMULATION 

3.  SIMULATE  TIME  SERIES  FROM  AMPLITUDES  (USES  OUTPUT  CRQM 

STEP  Hi.  H 2.  OR  #4..) 

8.  HELP 


**********  4 


(PLEASE  ENTER  YOUR  CHOICE  AND  KEY  RETURN) 


DO  YOU  WANT  TO  SAVE  THE  OUTPUT  "ROM  LAST  OPTION 
STEP  On  A  SEPARATE  USER  SPECIFIED  FILE  ?  Y  OR  N 


Pro o rammed  S'CP 


Time  series  are  stored  in  printer 
listing  UST45.DAT 


-:SIM7 


CONSOLE  UST 
OPTION  NO.  7 


**************************************************** 

PROGRAM  SIMULATES  WAVE  PROPERTIES  BY  FREQUENCY  DOMAIN  METHOD 
FOR  EITHER  UNCONDITIONAL  OR  CONDITIONAL  SIMULATIONS. 
(WRITTEN  BY  LEON  BQRGMAN,  LARAMIE,  WYOMING) 

PLEASE  KEY  RETURN  TO  CONTINUE 

********************************************************** 


DO  YOU  WISH  TO  READ  OUTPU"  FROM  PREVIOUS  RUNS 
AS  INPUT  -OR  THIS  RUN?  (Y  =  YES.  N  =  NO 

Y 

WHAT  IS  THE  FILE  NAME  FOR  THIS  INPUT? 

J :  OU'r 4.  DAT  -  INPUT 

OPTIONS  ACTIVE  IN  THIS  SUBPROGRAM: 

0.  EXIT  PROGRAM 

7.  PRODUCE  A  .REDUCED  SET  OF  AMPLITUDES  FOR  A  LOW  RESClUT I C  '■) 
REPRESENTIQN  THE  SEA  SURFACE  AND  KINEMATICS. 

(USES  OUTPUT  FROM  S"E^'S  *1  OR  S^EP  *3.  ) 
a.  HELP 

************************************************** 

(PLEASE  KEY  YOUR  CHOICE  AND  RETURN) 

7  - OPTION  NO. 

PLEASE  ENTER  THE  LOWEST  AMPLITUDE  OF  INTE.RES- 

0.  1 

QWKSRT 

DO  YOU  WANT  TO  SEE  A  COMPARISON  OF  THE  FULL 
RESOLUTION  AND  THE  LOW  RESOLUTION  WAVE  TIME  HISTORIES 
AT  X-0.  Y»0  ?  <Y  OR  N) 

Y 

ENTER  TIME  STEP-  FOR  START  OF  LIST  AND  RETURN 
1 

ENTER  TIME  STEP  AT  TERMINATION  uF  LIST  AND  RETURN 
300 

ENTER  0  IF  GRAPHICAL  OUTPUT  IS  DESIRED 
KEY  1  IF  NUMERICAL  LIST  OF  TIME  SERIES  IS  WANTED 
0 

OPTIONS  ACTIVE  IN  THIS  SUBPROGRAM: 

0.  EXIT  PROGRAM 

7.  PRODUCE  A  REDUCED  SET  OF  AMPLITUDES  FOR  A  LOW  RESOLUTION 
REPRESENTIQN  THE  SEA  SURFACE  PND  KINEMATICS. 

(USES  OUTPUT  FROM  STEPS  #1  OR  S"^-P  #3.  ) 
a.  HELP 

e************************************************* 

(PLEASE  KEY  YOUR  CHOICE  AND  RETURN) 

^r,  vr*"  .ion'*’’  ■  »cr  n»  itoi  it  er?niv*  •  •*>0*r*  r* -•  -r - m 


* 


Programmed  STOP 
F> 


List  of  low  resolution  amplitudes  ie 
printed  on  line  printer  as  shown  in  UST7.DAT. 


DISTRIBUTION  LIST 
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